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We investigate a population-imbalanced two-species fermionic system where the 
resonantly-paired fermions combine to form bosonic molecules via Feshbach inter- 
action. The natural dynamics of the system is studied and it is shown that the 
oscillation of the condensate fraction is periodic or quasi periodic, depending on 
the value of Feshbach coupling. We describe how a time dependent magnetic field 
can be used to study the natural frequencies and thus explore the momentum space 
structure of the population imbalanced system. 
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I. INTRODUCTION 

Ultracold fermionic atoms have attracted a lot of attention in the last decade. The 
fermions here can exhibit BCS pairing or can form molecules and undergo Bose Einstein 
condensation (BEC). The crossover from BCS to BEC can be controlled by a magnetic field 
capable of providing Feshbach resonance lMl7|. An interesting variation in this situation is 
the introduction of an imbalance in the fermion number. Pairing now has to work around 
the fact that not all fermions of type A(one of the species) have a fermion of type B to 
pair with. This is equivalent to considering superconductivity in the presence of an external 
magnetic field which creates an imbalance of sp in up and spin down states. The consequence 
can be exo tic p airing states like FFLO 18], ll9], Sarma (breached pair) phase 20H23| , phase 



separation 



16 



24j . These are, however, notoriously difficult to detect experimentally. 
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Dynamics is a very important attribute of ultracold systems and has been studied exten- 
sively. R.A. Barankov et al. discussed the collective nonlinear evolution of BCS state when 
the pairing interaction is turned on instantaneously 25|. V. A Andreev et al. showed that 



when the position of the Feshbach resonance is changed abruptly, the superfluid undergoes 



a coherent BEC-to-BCS oscillation 26]. M. H. Szyman'ska et al. |27j addressed the situation 
when the Feshbach magnetic field jumps suddenly, and studied the short time dynamics that 
followed - damped oscillations with an amplitude depending on initial conditions. They also 
argued that atom-molecule oscillations are negligible for all practical purposes. All these 
authors concentrated on following the system after a rapid change. This entails the study 
of the nonlinear term in the evolution equation. 

In this paper, we aim at studying the dynamical properties of a population-imbalanced 
Fermi gas but without introducing a sudden change in Feshbach magnetic field. We focus 
on the linear dynamics of the system when a natural fluctuation of BEC condensate occurs 
along the BCS-BEC crossover path, and the system tries to relax back to the steady state. 
In the process we find that it shows periodic or quasiperiodic oscillations. The strength of 
the Feshbach term determines whether the oscillation will be periodic or quasi periodic. The 
frequencies of motion are found to be sensitive to the pairing in the momentum space. We 
propose that a time dependent Feshbach coupling can be used to probe the pairing structures 
in momentum space. This can provide another handle on exploring experimentally the exotic 
pairing states. 

I. MODEL HAMILTON! AN 

Here we start with a two-species fermionic system. In addition to the fermion-fermion 
interaction ( denoted by gi), there is an additional interaction (g 2 ) of the Feshbach variety 
which couples two fermions to form a bosonic molecule . Our model resembles the one used 



in 
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28l | . ijj p . and ifi- P ^ denote the fields for the two fermions, and (f> represents the 



Bosonic field. 

Our Hamiltonian is : 
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p p v 

(1) 

Next we calculate the commutation relations in order to arrive at the equations of mo- 
tion. This is in line with the Ehrenfest Theorem, which relates the time derivative of the 
expectation value for a quantum mechanical operator to the commutator of that operator 
with the Hamiltonian of the system. 

Heisenberg operators are introduced as t^H^ipt) = e lHt / h i(jp^e~ lHt / h and ^n\{pt) = 
giHt/h^J e -iHt/n _ fjkg e q Ua tions of motion are : 



(2) 



Let us now define the pair wavefunction. 

P = (^TV^> (3) 

It turns out that 

80 

ih-^ = 2e p (i(, H ^ Hl ) + g(iPM(N p - 1) + g 2 <P(N p - 1) (4) 

Here N p — N\ + N 2 , i.e, the total population corresponding to a particular momentum p. 
Now, (ipH^ipHi) represents the expectation value of the pair wave function in the Heisenberg 
picture while (V'tV'j.) * s the same quantity in the Schroedinger's picture. Since these are 
expectation values, they should not depend on the choice of the representation. So we can 
replace (i^h^hi) by O p as well. Moreover, when both the pairing states are occupied, (e.g, 
in region of BCS pairing), N p = 2 Therefore, 

dO 

ih-^- = {2e p + g)O p + g 2 <t> (5) 



As for the evolution of 

dt 



ih ~ftF = 92^2 O p + e b c/) (6) 



p 
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A quic k comparison shows that these are linearized versions of the evolution equations 



in 



26 



27] , which is the relevant part for our calculation. The non linear terms are essential 



for dynamics following a sudden quench. 

II. SYSTEM IN EQUILIBRIUM: 

The first test of the new theory is that, whether we can reproduce previously established 
results as special cases of the new formalism. So we try to go back to the static case using 
these dynamical equations. 

We define a single particle Green's function in the momentum space. 
g(pt,p't f ) = -([^(pt^H^ip'f)}) and P{pt,p't') = -([^H t \pt)^H t \p't')}) and find 
that 

ihj t g(pt,p't') = ih(p - p')5(t - H) + e p g(pt,p>t>) + g eff {^ n )^(pt,p't') (7) 



A direct analogy with the standard BCS 29[ leads to the relation 



9eff X)<^MW> = sE^tW + 92<f> (8) 

V V 

In the static case, let the condensate order parameter be a constant, i.e, — = 

at 

g 2 

From Equations (EI) and (jHJ), g e ff — 9 -If the fermion-fermion four point interaction is 

an attractive one (which is indeed the case for a BCS superfluid), then taking The effective 
interection g e ff in the attractive sense and writing g = — \gi\, we arrive at 

9e f f = gi + g f (9) 
tb 



This expression matches with the one obtained using diagrammatic methods H[ and 



variational technique [28] So, here we are able to reproduce the same results without having 
to use a trial form of the ground state of the system. The Ehrenfest theorem and the 
commutation relations have directly landed us here, and this treatment is certainly more 
general in nature. 
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III. SYSTEM OUT OF EQUILIBRIUM : FREQUENCIES OF OSCILLATION 

Let, O p be the fluctuation in O p , and <fi be the fluctuation in 0, Therefore, 

ih-^ L = (2e p + g)O p + g 2( j) (10) 

ik— = g 2 s ^6 p + e b (j) (11) 
k 

g 2 (j>(uj) 



We take the respective Fourier transforms and find that O p (u) = 

a p + nu 

So far, the treatment has been general, and the population balance or imbalance has 
never been taken into account. Now we come to the specific situation where there are two 
species of Fermions, or may be, two hyperfine states of the same atom (for simplicity, we 
have used the t and 4- to denote them): but the one has a larger population than the other. 
In the momentum picture, this corresponds to a particular geometry. The natural choice is 
that of a two-tier structure, one consisting of the paired superfluid, and another comprising 



of a normal Fermi liquid made of the remaining unpaired fermions. It is shown that 28] a 
population-unbalanced Breached Pair state is stable only when the core is normal and the 
outer shell, superfluid. 




Superfluid 



FIG. 1: Stable Structure for Population Imbalanced 'Breached Pair' State 

In the above figure, the shell structure in momentum space is depicted. The unpaired 
majority fermions stay in the core region, which is a normal fluid. The pair superfluid forms 
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the outer shell, from momenta pi to p 2} where pi and P2 form the momenta boundaries 
between which the fermions form BCS-like pairs. 

Experimental results obtained by the Ketterle group at MIt[s^], too, confirms this struc- 
ture. In a strongly interacting Fermi gas with imbalanced populations, they observed a 
superfluid region surrounded by a normal gas in the form of a shell structure in the coordi- 
nate space. Now, the structure we described above (normal core and superfluid shell), when 
mapped into real space via a Fourier Transformation shows a high density of superfluid in 
the centre ( i.e, the picture obtained by Shin et al.) provided the population imbalance is 
not too high. Our Fourier space calculation shown that this holds till p 2 > 1.26pi. Now, 21 
is the measure of the population imbalance of the system, and indeed from the experimental 
results obtained by Shin et al., this core and shell structure survives upto a population im- 
balance of 75 percent (i.e, when p 2 > 1.05pi). So our calculation matches with their findings 
within an error margin of 20 percent. 

Thus, O p (u) has to be summed over the superfluid region, i.e, over all p from p\ to p 2 . 

?^ = -**f(OT (12) 

Let pf be the Fermi momentum and M the mass of the atoms. Therefore a p + uh = ap + b, 
where a = b = —-jf + g + ojH. Integrating Eq.( IT2|) by parts, we obtain ^2 p O p (uj) = 
-g 2 <f>(uj)(f(p 2 ) - /( Pl )) where f(a) = - §ln(p + J) - ?(p + J) + fln(p + J)] 

Putting back in Eq fllll) . we obtain 

<f>{u})[e b + tuj-gl{f(jp 2 )-f{ Pl ))] = Q 

Which means, 4>(u) is zero if [eb + hu — g\j(j>2) — /(pi))] 7^ 0. Therefore, in the expansion 
of (j)(t), only those (fi(u)s will survive for which 

e b + f Uj j-g 2 2 (f(p 2 )-f(p 1 )) = (13) 

Therefore, <p(t) = (pie iuJlt + <p 2 e iuJ2t + ... 

Here u\, ui 2 .. are the solutions of equation (TI31) . 

To find out wheter there exists real values for u, we take resort to graphical solutions. We 
are only interested in real co, because that would give us solutions in the form of 4>{t) = (f}e tult , 
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which denotes oscillation. If, on the other hand, we get imaginery solutions for co, then the 
solutions are of the form <p(t) = 0e w * or 0(£) = 0e _iW *. The first one signifies an exponential 
growth in <p and is, therefore, unphysical. The second one marks exponential decay, and its 
effect should be negligible as time increases. So we look out for real solutions only. 

From equation f|T3|) . the roots of cu are given by u = fl(uj), where fl(ui) = g2i.fi.P2) — 
/(Pi)) — e b- So we plot u along X-axis and both fl{ui) and to along Y-axis. The blue lines 
correspond to and the red lines mark u. If the two lines cut at any point, we get real 

solutions of the equation at that point. 

We scale all energies by the Fermi energy Ep, and all momenta by the Fermi momentum 
Pf of the majority species. Therefore, in this convention, mass of each partcle gets fixed at 
0.5. h is taken as 1. 

When the Imbalance is Fixed and $2 is Varied: The BCS pairing takes place near 

the Fermi surface, within a cut-off region. For standard superconductors, this cut-off is ftwn, 

4 \ 

ujd being the Debye frequency. For Ultracold atoms, the cut-off is — Ep=0.5AlEp 31]. 
Since we have scaled all energies by Ep, Fermi level corresponds to 1, and the lower cut-off 
is at 1 — 0.541 = 0.459. Thus we have to choose pi to lie between 0.459 and 1. We take P2 
to be 1, the Fermi momentum. 

We take p\ = 0.5, P2 = 1, eb at 0.03. From figure 5.2, we see that when g 2 is 0.1, there 
is only one point where /l(w) and ui cut one another, i.e, only one real solution for u>. A 
single solution exists for #2 = 1 as well. As this coupling is increased slightly, at c/2 = 1-7, 
there appears two such points, i.e, two real us. Then, as ^2 obtains higher values, there are 
always two real solutions for u . 

So we can call g 2 = 1.7 a critical coupling. If g 2 is less than this coupling value, there is 
only one real frequency of oscillation in 4>{t). Beyond it, there are two frequencies. 

When (72 is Fixed and the population imbalance is Varied: Now P2 has been fixed 
at 1, 6b at 0.03. We vary pi, which measures the amount of imbalance, because it is up to 
momentum p\ that the majority fermions remain unpaired. 

In figure 5.3, we observe that when the population imbalance is very low (pi is 0.5, i.e, a 
value slightly higher than the lower cut-off for pairing), the double frequency region appears 
at g 2 = 1.7. As the imbalance is increased gradually, the value of the critical coupling 
decreases. 
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FIG. 2: variation of fl(oj) with oo and solutions for fl(co) = to for different values of g 2 

Thus, the fluctuation in <fi can undergo oscillation with 1) one frequency 2) two frequencies 
3) no frequency at all, depending on the value of g 2 . 

In Case 1, <j>(t) = <\) X e ibJ ^\ or, </>(t) = (p + ^e iu)lt 

In Case 2, <j>(t) = fae 1 ^ + (p 2 e iuJ2t , or, <p(t) = (j) + (f)ie iuilt + (f) 2 e iuJ2t 

In Case 3, = cf>o, because the fluctuation decays exponentially. 



IV. PROBING BY AN OSCILLATORY DRIVE 



Let us add a small oscillatory component to the magnetic field. So, instead of g 2 , the 
coupling is of the form g 2 (l + ee lQt ) where e << 1. We can make a perturbative expansion 



4>{t) = 4>° + 14>' 

0(t) = 0° + eO' 
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a) a) 

FIG. 3: variation of with u and critical coupling for different values of population imbalance. 

Where cjP and 0° are the values of <j>{t) and 0(t) when there is no oscillatory part in the 
coupling. Noting that cjP = 4>ie luJt , it follows that 

BO' 

ih— = a p O' + g 2 (j>' + gzfaeV**** (14) 
Taking Fourier Transforms as before, 



,// \ gl^uj^-KdiVL + uj x - u)(f(pi,p 2 ,u 1 ) + f(p 1 ,p 2 ,u)) 
<P M = 7- VT, \ (15) 

e b + fiw -gif(p,u)) 

Here f(pi,p 2 , w) = f(p 2 , oj) — f{p\, u>). Now, 4>'(u) is non-zero only when co = u)\ + Q and 
its value at that particular frequency is 

,// x = g%M UJ )2n(f(pi,P2,Ui ±ty + f(pi,P2,Wi)) (m) 
9K ' e b + h(n + oo l )-glf(p 1 ,p 2 ,u; 1 + n)) 1 ' 

There is a resonance when the denominator becomes zero, i.e, e& + h(Q + u\) — 
g^f(pi,p 2 ,ui + Q) = But we know, e& + hw\ — g 2 f(pi,P2,^i) = 0. Subtracting, 

KM + gl{fiPl,P2,Ui) - f(pi,P2,Wl + fi)) = (17) 

Since Q is associated with the Feshbach resonance, by tuning the frequency of the time- 
dependent magnetic field, we can control Q. If, for a particular fl the above equation is 
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satisfied, then we have a sharp resonance in the fluctuation is 0. 

If, <p° = <pie luJlt + (f) 2 e lul2t , following the same treatment, (f)'(oj) is non-zero only when 
oj = {ui\ + VL) or (u>2 + fi), and its values at those frequencies are 

g%(f) l (u)27r(f(p 1 ,p 2 ,uj l + n) + f(p l ,p 2 ,uj l )) ^ gl4> 2 (u)2n(f(p ) u 2 + g) + f{pi,P2,u 2 )) 
tb + hiSl + ux) - glf(p, Ui + Q)) e b + h(tt + u 2 ) - g 2 f(pi,p 2 ,u 2 + fi)) 

respectively. 

Therefore, resonance will take place when f2 satisfies either of the equations: 

fXl + 92(f(Pl,P2,Ul) - f(Pl,P2,Ul + ^)) = 

(18) 

^ + ^2(/(Pl>P2,W 2 ) - f(p!,p 2 ,UJ 2 + U))=0 

It is obvious that f2 = is a trivial solution of the equations. Again we take resort to 
graphical solutions to find out whether non-zero real solutions exist for Q. 
Let g$(f(pi,p2, ,wi) - /(pi,P2,w + fi)) = /'(^) 

So we plot f2 along X-axis and both /'(fi) and along Y-axis. The blue lines correspond 
to f'(Q) and the red lines mark Q. If the two lines cut at any point, we get real solutions of 
the equation at that point. 

Experimentally one can detect the oscillation frequencies. If an oscillatory component is 
added to the Feshbach coupling, and its frequency is tuned, then at a particular frequency 
there is a sharp resonance in the condensate fraction. We propose an algorithm for deter- 
mining the imbalance: chose a p\ (p 2 is set to 1 because at zero temperature it is the Fermi 
momentum of the majority species, by which all momenta are scaled)so that the numerically 
computed oscillation frequencies match with those obtained in the experiment. Then use 
Those frequencies and the external resonant frequency to calculate p\. If this doesn't match 
with the initial pi, use the new values as input for the first step. Iterating the process a 
number of times, one should arrive at the actual value of the breaching point, or the point 
which marks the separation in the momentum space. 
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FIG. 4: variation of /'(^) with and solutions for f'(£l) = £1 for different values of gi and 
corresponding u>. 

V. SUMMARY AND DISCUSSION 

Here we have studied the natural dynamics of a population-imbalanced fermionic system 
capable of making a BCS-BEC crossover. We have shown that the condensate fraction shows 
a periodic or a quasi periodic oscillation, depending on the value of the Feshbach coupling. 
There is a critical coupling below which the oscillation is always periodic, and its value 
depends on the population imbalance. 

We have shown that if there is an oscillatory component in the Feshbach term, one can 
achieve a sharp resonance in the condensate fraction by tuning the frequency of the external 
magnetic field. The breaching momentum can be calculated from this frequency value. 
Thus it proves to be an indirect method of experimentally determining the momentumspace 
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structure of the imbalanced Fermi system. 

This treatment can be extended to detect more complicated structures in the momentum 
space, too. The FFLO state, which arises due do finite momentum pairing, should also get 
reflected in a similar dynamical study. The finite pairing momentum is then embedded in the 
Bosonic field. If we use a single wave number q in operator <fi in our original Hamiltonian 
to represent the finite momentum, the oscillatory dynamics and the resonant frequencies 
should contain information regarding q, and thus portray both the momentum structure 
and the nature of pairing. 
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